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THE STABILITY OF ECCENTRICALLY 
STIFFENED CIRCULAR CYLINDERS 

VOLUME VI 

INTERACTION BEHAVIOR 
By 

G. W. Smith* E. E. Spier, 

S. N. Dharmarajan, and L. S. Fossum 
General Dynamics Convair Division 
San Diego* California 


ABSTRACT 

This is the last of six volumes, each bearing the same report number, 
but dealing with separate problem areas concerning the stability of eccen- 
trically stiffened circular cylinders. The complete set of documents was 
prepared under NASA Contract NAS8-11181. This particular volume includes a 
digital computer program which determines critical combinations of axial 
load (tension or compression) and radial pressure (bursting or crushing) . 

The theoretical basis for this program is a solution developed by the NASA 
Langley Research Center. No consideration is given to the influences of pre- 
buckling distortions from the true cylindrical shape. Under the sponsorship 
of the NASA Marshall Space Flight Center, the General Dynamics Convair 
Division has extended the foregoing theoretical solution to incorporate the 
influences of external running shear acting either alone or in combination 
with axial load and/or radial pressure. The extended solution is also 
presented here. In addition, some discussion is given of the case where 
eccentrically stiffened cylinders are subjected to overall bending acting 
either alone or in combination with axial load. 
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Symbol 



A 

r 


AAA 
11’ 12’ 13’ 

A 22 ,A 23’ A 33 


a 


a 


n 


B 

m 


b 

n 


C 

m 


c 

n 


D 

m 


DEFINITION OF SYMBOLS 
Definition 

Amplitude of buckling displacement component 
[see equations (2-12)]. 

Cross-sectional area of single ring (no 
cylindrical skin included). 

Cross-sectional area of single stringer 
(no cylindrical skin included). 


Parameters defined by equations (2-ll). 
Overall length of cylinder. 

Amplitude of buckling displacement component 
[see equations (2-18)]. 

Amplitude of buckling displacement component 
[see equations (2-12)]. 

Amplitude of buckling displacement component 
[see equations (2-18)]. 

Amplitude of buckling displacement component 
[see equations (2-12)]. 

Amplitude of buckling displacement component 
[see equations (2-18)]. 

Amplitude of buckling displacement component 
[see equations (2-12)]. 
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Symbol 



D 

y 


D 

xy 


d 

E 

E 


E 

r 


E 

s 


E , E 
x’ y 


F 

m 


G 

G(m) 

G 

r 

G 

s 

G 

xy 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Bending stiffnesses of basic cylindrical skin in 

longitudinal and circumferential directions, 

respectively. (For isotropic skins, 

D = D * Et 3 /12). 
x y 

Twisting stiffness of basic cylindrical skin 

3 

(For isotropic skins, D = Gt /6). 

xy 

Stringer spacing. 

Young's modulus. 

Amplitude of buckling displacement component 
[see equations (2-12)1. 

Young's modulus for ring. 

Young's modulus for stringer. 

Extensional stiffnesses of basic cylindrical skin 
in longitudinal and circumferential directions, 
respectively. (For isotropic skins, E = E = Et). 

x y 

Amplitude of buckling displacement component 

[see equations (2-12)1. 

Modulus of elasticity in shear. 

Parameter defined by equation (2-15). 

Ring modulus of elasticity in shear. 

Stringer modulus of elasticity in shear. 

In-surface shear stiffness of basic cylindrical 

skin (For isotropic skins, G = Gt). 
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Symbol 

I 

o 


I 


o 

s 


I 

r 


I 

s 


J 

r 


J 


s 


M ,M , M 
x’ y’ xy 


m 


N , N , N 
x’ y’ xy 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Moment of inertia of single ring (no cylindrical 
skin included) taken about middle surface of basic 
cylindrical skin. 

Moment of inertia of single stringer (no cylindrical 
skin included) taken about middle surface of basic 
cylindrical skin. 

Centroidal moment of inertia of single ring (no 
cylindrical skin included) 

Centroidal moment of inertia of single stringer 
(no cylindrical skin included). 

Torsional constant of single ring (no cylindrical 
skin included) . 

Torsional constant of single stringer (no cylindrical 
skin included). 

Ring spacing. 

Bending and twisting stress resultants of basic 
cylindrical skin. 

Number of axial half-waves in buckle pattern; 
summation index. 

Normal and in-surface shear stress resultants of 
basic cylindrical skin. 
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Symbol 

N 

x 





N 

*y 


N 

y 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Uniformly distributed applied longitudinal 
running load acting at the centroid of skin- 
stringer combination (positive for compression; 
negative for tension) . 

Peak value of non-uniformly distributed component 
of applied longitudinal running load acting at the 
centroid of skin-stringer combination [see equation 
(2-17)]; (positive for compression; negative for 
tension) . 

Value of uniformly distributed component of applied 
longitudinal running load acting at the centroid of 
skin-stringer combination [see equation (2-17)]; 

(positive for compression; negative for tension). 

See equation (2-17). 

Critical value of uniformly distributed longitudinal 
running load N^ when acting alone. 

Uniform running shear load due to applied overall 
external torque acting about axis of revolution. 

Uniformly distributed applied circumferential running 
load acting at the centroid of the skin-ring com- 
bination (positive for compression; negative for tension). 
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Symbol 


N 

y 


o 


n 


R 

R 


R 

y 


R 1 or a 2 


t 

u,v f w 


U, V, w 


x,y,z 


z 


r 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Critical value of uniformly distributed circum- 
ferential running load N when acting alone. 

y 

Number of circumferential full-waves in buckle 
pattern; summation index. 

Radius to middle surface of basic cylindrical skin. 

Ratio of applied axial loading (N^) to the critical 

value of axial loading when acting alone (N ) 

_ o 

Ratio of applied circumferential loading (N ) to the 

y 

critical value of circumferential loading when acting 

alone (N ) . 
y o 

Ratio of an applied load (or stress) to the critical 
value for that type of load (or stress) when acting 
alone. 

Thickness of basic cylindrical skin. 

Displacements in the x, y, and z directions, 
respectively, of a point in middle surface of basic 
cylindrical skin. 

Amplitudes of buckling displacements [see equations 
( 2 - 9 )]. 

Longi tudinal^ circumferential ? and radial directions, 
respectively. 

Distance from centroid of ring (no basic cylindrical skin 
included) to middle surface of basic cylindrical skin 
(positive for external rings; negative for internal rings). 
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Symbol 


MM 



^Axial 

r Bend 


6 . 
^pm 


c ,e , y 
x’ y 'xy 




V ■ 

P * l 1 
x y 





X 


c 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Distance from centroid of stringer (no basic 
cylindrical skin included) to middle surface of 
basic cylindrical skin (positive for external 
stringers; negative for internal stringers). 
Knock-down factor for circular cylinder subjected 
to pure axial load (see Volume V). 

Knock-down factor for circular cylinder subjected 
to pure bending (see Volume V). 

Quantity defined by equations (2-2l). 

Quantity defined by equations (2-16). 

Strains at middle surface of basic cylindrical skin. 
Longitudinal strain of stringer. 

Circumferential strain of ring. 

Poisson's ratio (For isotropic material). 

Poisson's ratios for bending of orthotropic skin in 
longitudinal and circumferential directions, 
respectively (For isotropic skins, 



Poisson's ratios for extension of orthotropic skin 
in longitudinal, and circumferential directions, 
respectively. (For isotropic skins, 

p = p =p # =p'=p). 
x y x y 

Total change (due to buckling displacements) in 
potential energy of loaded stiffened cylinder (Also 
used as conventional notation for the constant 
3.14 ). 

Change (due to buckling displacements) in strain 
energy of basic cylindrical skin. 
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Symbol 

*L 


it 


r 


o 

cr 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Change (due to buckling displacements) in 
potential energy of external loading. 

Change (due to buckling displacements) in 
strain energy of rings. 

Change (due to buckling displacements) in strain 
energy of stringers. 

Critical buckling stress 


NOTE: 


A subscript preceded by a comma denotes partial differentiation 
with respect to the subscript variable. For example, 


w, 


xy 


>2 

q w 

5x5y 


xv 
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SECTION 1 
INTRODUCTION 

For isotropic cylinders that are subjected to combined external 
loads, it is customary to represent the critical loading combinations 
by means of so-called interaction curves. Figure 1 shows the graphical 
format usually employed for this purpose. 



Figure 1 - Sample Interaction Curve 
for Isotropic Cylinders 

The quantity R^ is the ratio of an applied load or stress to the critical 
value for that type of loading when acting alone. The quantity R^ is 
similarly defined for a second type of loading. Curves of this type give 
a very clear picture as to the structural integrity of particular configura- 
tions. All computed points which fall within the area bounded by the inter- 
action curve and the coordinate axes correspond to stable structures. All 
points lying outside this region indicate that buckling will occur. Further- 
more, a measure of the margin of safety is given by the ratio of distances 
from the calculated point to the curve and the origin. In view of the 
history of satisfactory experience with this type of curve in the design 
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and analysis of isotropic cylinders, it is quite natural to consider the 
application of this concept to eccentrically stiffened cylinders that 
are subjected to combined loads. However, whereas the behavior of isotropic 
cylinders can be adequately covered by a single curve for each two-load 
combination, the findings reported in this volume show that a multiplicity 
of curves is required for each two-load combination applied to eccentrically 
stiffened circular cylinders. This complexity arises out of the greater 
number of independent variables required to describe the physical features 
of the latter configurations. To keep the number of such curves within 
reasonable bounds, one must resort to simplifications which result, of 
course, in a loss of rigor. However, sufficient accuracy can be retained 
to make such curves useful for initial sizing, rough checking, and the study 
of trends. For the purposes of a final design, more refined values would 
be required. In general, these can only be obtained as single-point solutions 
from digital computer programs. 
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SECTION 2 


EQUATIONS 


2.1 COMBINED AXIAL LOAD AN!) RADIAL PRESSURE 


The theoretical solution by Block, Curd, and Mikulas [ l] provides 
a suitable means for the stability analysis of eccentrically stiffened 
circular cylinders that are subjected to combined axial load (tension or 
compression) and radial pressure (bursting or crushing). This solution 
has been verified through the rederivation of reference 2. In general, the 
development by Block, et al. C ll was accomplished by the application of 
variational techniques to establish bifurcation points (see Volume I [3]) 
along an initially linear equilibrium path. This was achieved by first 
formulating expressions for the changes in strain energy due to the buckling 
displacements. For the basic cylindrical skin this increment was expressed 
as follows: 

2rcR a 

it=i( I (N e + N y + N c - M w 
c 2 J J xx xy ’ xy yy x,xx 

0 0 ( 2 . 1 ) 

+ 2M w - M w ) dx dy 
xy ,xy y ,yy 


where 



xy 

xy 



Stress resultants of basic cylindrical skin. 


c , e , Y = Strains at middle surface of basic cylindrical 
X y Xy skin. 

= Change (due to buckling displacements) in 
strain energy of basic cylindrical skin. 

= Longitudinal direction. 

= Circumferential direction. 

= Cylinder radius. 

= Overall length of cylinder. 

Numbers in brackets [ ] in the text denote ref erences listed in SECTION 6. 
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The manner in which Block, et al [l] formulate the stress resultants 
facilitates analysis where the basic cylindrical skins themselves have 
orthotropic properties. Preceding then to the longitudinal stiffeners, 
their change in strain energy was expressed as follows: 




dA 

s 



( 2-2 


where 


A 

s 



d 

G 

s 

J 

s 

w 


= Change (due to buckling displacements) in strain 
energy of longitudinal stiffeners. 

= Cross-sectional area of longitudinal stiffener (no 
cylindrical skin included) . 

= Longitudinal strain of longitudinal stiffener. 

= Stringer spacing. 

= Shear modulus of longitudinal stiffener. 

= Torsional constant of longitudinal stiffener. 

= Radial displacement. 


The change in strain energy of the circumferential stiffeners was found 
from 



where 

l = Ring spacing 

r = Subscript denoting ring (circumferential stiffener 
not including any cylindrical skin). 


(2-3) 


2-2 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC-DDG-67-006 


The change in the potential energy of the external loading was taken as 


n L 



+ 2N w. w, 
xy ’ x ’y 


+ N w, 

y 



dx dy 


(2-4) 


where 


*L 


Change (due to buckling displacements) in potential 
energy of external loading. 


= Applied running axial load. 

= Applied running hoop load. 

_y 

N = Applied running in-surface shear load. 

Although this expression includes the applied running shear load N , 

xy 

Block, et al. [ ll later set this term equal to zero so that their final 
equation applies only to the case of applied axial load and/or applied 
radial pressure. 

By making use of equations (2-1) through (2-4), the total potential 
energy of the system was expressed as follows: 


it = it + it + it +it T (2-5) 

c s r L 

The next step was to employ the principle of stationary potential energy 
(see Volume I [3]) to arrive at the following set of equilibrium equations 
in terms of the buckling displacements u, v, and w: 


E 


1 - p ’p 1 

x y 


u, 


xx 



- 2 W, 

XX S XXX f 


+ 


1 - p 'p ' 
X y 



0 


(2-6) 
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E 


V, \ E A 


w, 

1 V « 


r-"v'{r-( v ’yy + -/) + - TT Ky + ‘ v» 


yyyj 


1i 'E 
*—X 


1 - i» x 'ii * ’xy xy \ ’xy ’xx 


^ U ’xy + v, xx] 


* — , u, + G ( u, + v» | = 0 

•It ' ' W VV 1 ' w vv I 


(2-7) 


D Vo \ / p D 

* + ) w , J -J jls- + 2D + s-y + -aj + 


p D G J G J 

v ^ ~ r r 


L 1 - ' d / "’xxxx T l 1 - p„p„ T “"xy 1 - Pji„ d T l j^'xxyy 


x y 


x y 


x y 


D E r I o 

+ \ r-^yT * -T-^Km 


E , 

R(l - p 'p ' ) ( ^x u 

x y V 


, + v« + 

’x ’y 


*) 


E A 
s s — 


E A _ / 2w, 

r r — 1 


j z u, — . « | n 

d s xxx £ r \ R 


yy + 


V 


\ E A . x 

1 r r / w \ 

’yyy / + R* ( v ’ y + r j 


+ N w, + 2N w, + N w, =0 
x ’xx xy ’xy y ’yy 


(2-8) 


Block, et al. Cl] then obtained a solution to these equations by assuming 
boundary conditions of classical simple support and the following set of 
buckling displacement functions: 


where 


u 


v 


w 


u cos 


v sin 


w sin 


mux 

a 

mux 

a 

mux 

a 


cos 


sin 


cos 


SX 

R 


SX 

R 


SX 

R 


(2-9) 


m = Number of longitudinal half-waves in buckle" pattern, 
n = Number of circumferential full-waves in the buckle 
pattern. 
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By substituting equations (2-9) into the equilibrium equations, it can 
be observed that the existence of non-trivial buckling displacements 
requires that a certain determinant vanish. This condition reduces to 
the stability equation 




+ / A 12 A 25“ A 13 A 22 
\ A 11 A 22" A 12 2 



A 12 A 13~ A 11 A 23 \ 
A 11 A 22" A 12 / 



( 2 - 10 ) 


where the A. ,'s are functions of the material properties, the geometry 
of the cylinder, and the shape of the buckle pattern. In particular. these 

7 

quantities can be expressed as follows: 
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Note that the foregoing A _ 1 s are no * elastic constants defined in 

Volume I [31 and used elsewhere throughout this multiple-volume report. 

liquation (2-10) is the fundamental buckling criterion that emerges 

from the derivation of Block, et al. C 1 1 - It should be observed that this 

equation does nothing more than establish N and N values which are capable 

x y 

of maintaining the cylinder in deformed configurations which correspond to 
particular numbers of longitudinal half-waves m and circumferential full 
waves n . Calling upon the bifurcation concept discussed in Volume I -3], 
it follows that the critical buckling loads can be established by exploring 
the possible deformed equilibrium configurations for a minimum load condition. 
The digital computer program of SECTION 5 makes use of equation (2-10) in 
precisely this manner. Through the input, the analyst prescribes the ranges 
of m and n to be explored. The machine then computes the or N values 
corresponding to each included m-n combination and prints out the lowermost 
load encountered. Further refinement to this program could easily be 
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accomplished to eliminate the prior judgements required of the user in 
selecting appropriate ranges of in and n to be screened. 

Although the foregoing solution can be considered suitable for the 
purposes of final design, it is not intended that the reader interpret 
this to mean that equation (2-10) is perfectly rigorous. At best, it re- 
presents a state-of-the-art operational capability in a rapidly changing 
technology. For one thing, it should be noted that this solution is based 
on small-deflection theory. Hence, the theoretical results obtained from 
equation (2-10) must be reduced by appropriate knock-down factors (see 
Volume V L4])to account for the detrimental effects of initial imperfections. 
Furthermore, the subject solution is based upon the complete neglect of non- 
cylindrical pre-buckling distortions which can be extremely important where 
pressure differentials exist across the shell wall. End and ring restraint 
to Poisson' s- ratio hoop growth further contribute to this type of pre- 
buckling deformation. It should also be observed that the solution of Block, 
et al. C l] is primarily based upon monocoque shell theory and its application 
to discretely stiffened cylinders is achieved by the conventional smearing- 
out technique. For this purpose, the stiffnesses of the discrete component 
members are averaged over the entire shell surface to obtain an equivalent 
monocoque analysis model. In view of this, it might sometimes be wise to 
consider the introduction of modified component section properties to arrive 
at true equivalence. Influences from such sources as shear lag, stiffener 
attachment techniques, etc. might be reflected into the analysis by such 
modifications. Finally, it should be noted that the solution of Block et al. 
[l] is based upon Donnell-type simplifications which render the results 
invalid for non- axi symmetric buckle patterns having a small number of cir- 
cumferential waves. The rule-of-thumb guideline is offered here that this 
solution be considered inapplicable for cases where O < n < 2. 

2.2 COMBINED AXIAL LOAD, RADIAL PRESSURE, AND RUNNING SHEAR 

To analyze this case of combined loading, the theoretical solution of 
reference 2 may be used. This solution constitutes an extension to the 
theory of Block, et al. [ll which is discussed in SECTION 2.1 above. To 
accomplish this extension, use was made of the same equilibrium equations 
as were employed by Block, et al. These expressions are given above as 
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equations (2-6), (2-7), and (2-8). Then, unlike the derivation of reference 

1, the N term was retained throughout all of the subsequent mathematical 

xy 

operations. The presence of this additional term made it necessary to 
introduce displacement functions which differ from those expressed as 
equations (2-9). In particular, the following formulations were selected: 


u = sin 


m=l,3 


m 


cos 


mtcx 

a 




£ cos 
ra 


mux 

a 


ra=2, 4, 


v = cos 


njL 

R 


B sin 


m=l,3, 


ra 


nwx ny O 

+ sm -fr / 

a R , 

m=2, 4, — 


( 2 - 12 ) 


w ss sin 



m=l,3, — 


C 

m 


sin 


nntx 


+ cos 


!SL V 

R 2-r 


m=2, 4, — 


D 

m 


sin 


nwx 

a 


These equations satisfy the boundary conditions of classical simple 
support. By substituting equations (2-12) into the equilibrium equations 
and by applying the Galerkin method [5], one can obtain the following set 
of homogeneous equations [2]: 


G(m) C 6 

m pm 


4N -g-V 
xy R Lji 


mp 


m=l,3, — 


m=2, 4, — 


2 2 

p -m 


D = 0 (p = 1,3,5, — ) (2-13 


m=2, 4, — 


G(ra) D 6 + 4N -l-V 

m pm xy R JLj 


m=l,3 , — 


— m j ? 3 ~ C m = 0 (p = 2,4,6, ) (2-1^ 

p -m 
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where 


G(m) = A, 


33 


A 23 (A 11 A 23 " A 13 A 12 ) * A 13 (A 22 A 13 " A 23 A 12^ 
A 11 A 22 " A 12 




(2-15) 


and 


6 =1 if p=m 

pm 

5 = 0 if p#m 

pm 


(2-16) 


The A^'s of equation (2-15) are the same quantities as were previously 
defined by equations (2-ll). 

Equations (2-13) through (2-16) furnish a theoretical basis for com- 
puting the buckling loads of eccentrically stiffened circular cylinders 
which are subjected to combined axial load and/or radial pressure and/or 
running shear. It should be noted that these relatively compact equations 
actually represent an infinite system of homogeneous algebraic equations 
and, from a practical standpoint, one must rely upon a digital computer to 
obtain numerical results. A typical application might involve the com- 
putation of a critical N in the presence of given N and N loadings. The 

xy x_ y 

desired critical value would be the lowest eigenvalue N for the determinant 
of the coefficients for C and D^. In general, the larger the selected 
determinant, the better is the accuracy but the greater is the computational 
effort (computer machine-time) . 

No digital computer program is furnished here to implement the fore- 
going analysis method. It is therefore recommended that such a program be 
developed in the near future. Note however that, since the above results 
constitute an extension to the basic development of Block et al. Cl], all 
of the limitations cited in SECTION 2.1 apply here as well. 

2.3 COMBINED AXIAL LOAD AND OVERALL BENDING 

In order to properly understand the state-of-the-art relative to this- 
combined loading condition, it is helpful to first consider the current 
status for stiffened cylinders subjected only to pure bending. To begin with, 
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it is pointed out that some disagreement exists among the recently 
published documents concerning this problem. On the one hand, the 
theoretical results of reference 6 indicate that the critical stresses 
for orthotropic cylinders under pure bending are equal to the critical 
values for pure axial compression. This is contrary to the findings of 
Hedgepeth and Hall [7] for a particular corrugated cylinder having in- 
ternal rings.. For pure bending, these theoretical results show a buckling 
stress which is approximately 1.23 times the buckling stress under uniform 
axial loading. Still another theoretical development was recently published 
by Block C 8] which similarly shows differences between the critical stresses 
for the two subject loading conditions. The theory developed by Block C 8] 
is an extension to the basic approach of reference 1. The same governing 
differential equations of equilibrium were employed. These expressions are 
given above as equations (2-6), (2-7), and (2-8). The values N and N were 

_ *y y 

each set equal to zero and was expressed in a form equivalent to 


N ) = N + N cos D 

x / x x. R 

c+b c b 


(2-17) 


This permits consideration of overall bending acting either alone or in 

combination with uniform axial loading N . The quantity N is the peak 

c x b 

running compressive loading due to the overall bending moment. In the 
solution of reference 8, the following expressions were used to describe 
the buckling displacements u, v, and w: 


u = cos 


mrtx 

a 


So 


C 

n 


cos 


tSL 

R 


v 


CD 


sin 



b 

n 


sin 



(2-18) 
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By substituting these functions into the equilibrium equations and 
applying the Galerkin method Cs], Block obtained the following final 
expressions: 



where , 


F 

n 



+ A 


23 


A 12 A 13 ~ A 11 A 23 
A 11 A 22 “ A 12 


+ A 


13 


A 12 A 23 " A 13 A 22 \ 
A 11 A 22 " A 12 2 /. 


(2-19) 


( 2 - 20 ) 


and 


n= 0, 1 , 2,3, » • • 


6 . =1 if ( j = n) 
jn 

6 . =0 if (j qfc n) 
jn 


( 2 - 21 ) 


The A _ ' s of equation (2-20) are the same quantities as were previously 
defined by equations (2-11). 

Just as in the combined loading analysis cited in SECTION 2.2 above, the 
relatively compact equations (2-19) through (2-21 ) actually represent an 
infinite system of homogeneous algebraic equations. Here again, from a 
practical standpoint, one must rely upon a digital computer to establish 
numerical values for the critical loadings. As before, the size of the 
selected system of equations will determine the accuracy obtained and the 
computational effort (computer machine-time) required. 

The basic nature of the numerical techniques involved in the application 
of equations (2-13) through (2-16) is the same as that associated with 
equations (2-19) through (2-21 ). To understand the fundamental approach, 
the latter set of equations will be further discussed. In this connection, 
it should be recalled that a set of homogeneous algebraic equations can 
have a nontrivial solution if, and only if, the determinant of the coefficients 
equals zero. Hence, in the case of equations (2-19) through (2-21 ), one 
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might first assemble the determinant of the coefficients for the a^'s. 

This determinant could then be expanded and set equal to zero to obtain 

a polynomial equation in which, for given m and n values, either N 

or N is the only unknown, F or any single m-n combination, the 
X b 

lowermost root of this equation could then be computed. Such computations 
might be repeated for all possible m-n combinations. The critical loading 
would be the lowermost value encountered in this screening type of operation. 
Although this procedure is relatively simple in concept, it is often quite 
difficult to actually determine the lowermost roots of the polynomial 
equations, particularly when the size of the selected determinant is large. 
Therefore, mathematicians have devised a number of numerical schemes to 
extract these roots (eigenvalues) from the determinantal equation. Most 
computer laboratories have on-the-shelf sub-routines to perform such 
operations and they will not be discussed any further here. 

Following his derivation of equations (2-19) through (2-21 ), Block [8] 
then computed some sample results for three contemporary types of eccentricallj 
stiffened circular cylinders: ring-stiffened corrugated cylinders, ring-and- 
stringer-stif fened cylinders, and longitudinally (stringer) stiffened 
cylinders. Like the results given in reference 7, Block C 8] obtained 
different critical stresses for the separate cases of pure axial load and 
pure bending moment. The ratios of 

Pure Bending 

Pure Compression 0 cr 

ranged from 1.013 up to 1.397. 

To place the contents of this particular section in a proper perspective, 
it might be noted here that the historical development of isotropic 
cylinder theory likewise involved some controversy as to the relative 
theoretical strengths for the separate cases of pure axial load and pure 
bending moment. An analysis presented in references 9 and 10 indicated that, 
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for an assumed buckle waveform, the critical stress for isotropic cylinders 
under pure bending was 1.3 times the critical stress for pure compression. 

This calculation was cited by Timoshenko [ll] without a qualifying statement 
as to the assumed buckle wavelength. Because test data seemed to substantiate 
the existence of such an increase, the 1.3 factor was used for decades as 
a generally applicable value. However, the small-deflection analysis of 
reference 12 recently revealed that the ratio of critical bending and 
compression stresses can vary widely with longitudinal wavelengths and 
that properly minimized results show that the correct ratio is essentially 
equal to unity. The apparent increase of bending strength over compressive 
strength indicated by isotropic test data can only be explained by a 
consideration of the sensitivity to initial imperfections. For axially 
compressed isotropic cylinders these defects result in a severe reduction 
of actual strengths below classical theoretical values. Since, under pure 
bending, only a small portion of the cylinder’s circumference experiences 
the peak stresses which initiate buckling, there is a statistical influence 
from the probability for defects to exist within this restricted region of 
the overall shell wall. As a result, it is reasonable to expect that, under 
pure bending, actual reductions below classical theoretical values will not 
be as severe as those for pure axial compression. 

In view of all the factors cited in this section concerning the inter- 
action behavior of stiffened cylinders and isotropic cylinders, it appears 
that, to facilitate the practical design and analysis of stiffened con- 
figurations, one might choose between the following two alternatives for 
cases of overall bending acting either alone or in combination with axial 
load: 

(a) Assume that the theoretical critical stresses are the same 
for the separate cases of pure bending and pure axial 
compression. The only differences between allowable levels 
for these two cases would then result from the application 
of knock-down factors (see Volume V [ 4 ]) which recognize 
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differences in the probabilities for initial imperfections 
to coincide with peak stress locations. The related inter- 
action curve would simply be taken as the straight line 
which intercepts the horizontal and vertical coordinate 
axes at values of unity. 

(b) Obtain theoretical critical loadings from digital computer 
solutions based on the formulation of Block [8]. As noted 
earlier, this solution gives different critical stresses 
for the separate cases of pure bending and pure axial com- 
pression. Here again, knock-down factors would be applied 
using the criteria of Volume V [ 4 ]. These factors would 
further accentuate the differences between pure bending and 
pure compression. To be structurally sound, a given con- 
figuration must be capable of supporting the loading com- 
bination 

/ Design N Design N \ 

i, M 

\ ^ Axial ^Bend / 

where 


^Axial = ^ noc k”down factor for circular 

cylinder subjected to pure axial 
load (see Volume V [ 4 ]). 

^Bend = Knock-down factor for circular 
cylinder subjected to pure 
bending (see Volume V E 4] ) . 

The recommendation of this volume is that alternative (b) be followed. 
However, no digital computer program is furnished here to implement this 
recommendation. Furthermore, note that, since the subject formulation of 
Block [8] constitutes an extension to the basic approach of reference 1, 
most of the limitations cited in SECTION 2.1 apply here as well. Only those 
limitations related solely to the presence of pressure differentials are 
inapplicable to the case under discussion. 


2-14 


GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG- 6 7-006 


SECTION 3 
PARAMETRIC STUDIES 


In order to observe some of the trends of interaction behavior for 
eccentrically stiffened circular cylinders, several studies were conducted 
for the case of combined axial load and radial pressure. Both of the 
following combinations were considered: 

(a) Axial compression plus radial crushing pressure. 

(b) Axial compression plus radial bursting pressure. 

These investigations were performed by using the digital computer program 
of SECTION 5. As noted earlier, this program is based on the theoretical 
solution of Block, et al. Cl]. It should be recalled that this solution 
completely neglects the effects of non-cylindrical pre-buckling deformations 
such as arise out of the presence of pressure differentials and discrete 
stiffeners. Localized restraint to Poisson's ratio hoop growth likewise con- 
tributes to these discontinuity-type deformations. These two influences 
are likely to be quite important for all stiffened configurations except those 
where the stiffeners are very closely spaced. 

The configurations studied here may be described in terms of the various 
input values to the digital computer program of SECTION 5. In particular, 
these are as follows: 


Configuration 1; 

R = Cylinder Radius = 38.6 in. 

a = Overall Length = 72 in. 

d = Stringer Spacing = 2.48 in. 

I - Ring Spacing = 6.00 in. 

E = 1.5x10® lbs/in. 

x 

D = 250 lb-in. 

x 

G = 2x10 psi 

xy 

p ' = 0.25 

x 

p = 0.30 

X 

E = 30x10® psi 

S 

G = 12x10 psi. 

s 

A = 0.020 sq. in. 

s n 


E = 2x10® lbs/in. 

y 

D = 300 lb-in. 

y 

D = 182.25 lb-in. 

xy 

u ' = 0.35 

y 

p =0.40 

E^ = 25x10® psi. 

G = 10x10® psi. 

r 

A = 0.040 sq. in. 
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I = 0.005 in. 4 

o 

s 

J = 0.004 in. 4 

s 

z =0 

s 

Configuration 2: 

Same as Configuration 1 except that 
f = Ring Spacing =0.5 in. 

Configuration 5: 

Same as Configuration 1 except that 
i = Ring Spacing = 72 in. 



0.010 in. 4 

0.006 in. 4 
0 


Note that the eccentricity values for the above configurations are all 
taken equal to zero which means that the stiffener centroids are located 
in the middle surface of the basic cylindrical skin. To provide some insight 
into the influences from non-zero eccentricities, the following configurations 
were also included in the study: 


Configuration 4: 

Same as Configuration 1 except that 

z = 0.50 in. and z = 0.75 in. 
s r 


Configuration 5: 

Same as Configuration 3 except that 

z = 0.50 in. and z = 0.75 in. 
s r 


Various critical combinations of applied axial compression (N^) and 
applied radial compression (N ) were determined for each of the preceding 

y 

configurations so that interaction curves could be plotted. The results 
are tabulated in Table I and are plotted in Figure 2 where the appropriate 
configuration numbers are shown in parentheses. The quantity R x is the 
ratio of applied axial loading (N ) to the critical value of axial loading 


when acting alone 


cumferential loading 


K ) a "' 


and the quantity R is the ratio of applied cir- 

_ y 

(N ) to the critical value of circumferential loading 

y 
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when acting alone ^ j. Observe that the interaction relationships 

cannot be expressed by a single curve. Even in the absence of eccen- 
tricities (configurations 1, 2, and 3), variations in the basic geometry 
led to different curves. The introduction of non-zero eccentricities 
resulted in further complication in that still a greater number of curves 
then emerged. 

To study the ’combination of axial compression acting along with radial 
bursting pressure, solutions were obtained for the same five configurations 
as cited above. These results are tfabulated in Table II and are plotted in 
Figure 3 where, once again, the appropriate configuration designations are 
shown in parentheses. The circumferential tensile loading due to internal 
pressure is shown nondimensionally in terms of the circumferential critical 
compressive loading j . Here too, it can be seen that variations in the 

basic geometry (including' eccentricity values) lead to different interaction 
curves. 
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TABLE I 


Calculated Data for Interaction 
Example Configurations - Axial 
Compression and External Radial Pressure 


CONFIGU- 

RATION 

N 

y 

— 

N 

X 

m 

n 

N v 

N 

x 

N 

y o 

N 

X 

o 

1 

0 

7868 

4 

7 

0 

1 


309 

7452 

3 

7 

.25 

.947 


567 

5901 

1 

5 

.459 

.75 


619 

5450 

1 

5 

.50 

.693 


791 

3934 

1 

5 

.639 

.50 


928 

2725 

1 

5 

.75 

.346 


1014 

1967 

1 

5 

.820 

.25 


1237 

0 

1 

5 

1 

0 

2 

0 

15157 

5 

5 

0 

1 


1761 

14536 

5 

5 

.25 

.959 


3521 

13915 

5 

5 

.50 

.918 


4983 

11368 

1 

3 

.708 

.75 


5282 

9931 

1 

4 

.75 

.655 


5699 

7579 

1 

4 

.809 

.50 


6371 

3789 

1 

4 

.905 

.25 


7042 

0 

1 

4 

1 

0 

3 

0 

4089 

1 

6 

0, 

1 


62 

3282 

1 

7 

.25 

.803 


74 

3066 

1 

7 

.301 

.75 


124 

2215 

1 

7 

.50 

.542 


185 

1148 

1 

7 

.75 

.281 


193 

1022 

1 

7 

.780 

.25 


219 

560 

1 

7 

.888 

.137 


247 

0 

1 

8 

1 

0 
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TABLE I 
(Continued) 


Calculated Data for 
Interaction Example Configurations - 
Axial Compression and External 
Radial Pressure 


CONFIGU- 

RATION 

N 

y 

N 

X 

m 

n 

N 

_2L 

N 

y o 

N 

x_ 

N 

X 

o 

4 

0 

8347 

3 

7 

0 

1 


266 

7396 

1 

5 

.25 

.886 


395 

6260 

1 

5 

.371 

.75 


533 

5049 

1 

5 

.50 

. 605 


632 

4174 

1 

5 

.593 

.50 


799 

2701 

1 

5 

.75 

.324 


868 

2086 

1 

5 

.816 

.25 


1065 

0 

1 

6 

1 

0 

5 

0 

4301 

1 

6 

0 

1 


62 

3405 

1 

7 

.25 

.792 


125 

2326 

1 

7 

.50 

.541 


187 

1247 

1 

7 

.75 

.290 


225 

564 

1 

8 

.90 

.131 


250 

0 

1 

8 

1 

0 
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y 




.8 1.0 


Figure 2 - Stability Interaction Results for 
Combined Axial Compression and 
External Radial Pressure 
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TABLE II 


Calculated Data for Interaction 
Example Configurations - Axial 
Compression and Internal Radial Pressure 


CONFIGU- 

RATION 

N 

y 

N 

X 

m 

n 

N 

N 

y 0 

N 

X 

N 

X 

° 

1 

123 7 

0 

1 

5 

i 

0 


0 

7868 

4 

7 

0 

1 


-309 

8202 

4 

7 

-.25 

1.042 


-518 

8535 

4 

7 

-.50 

1.085 


-1237 

9189 

4 

6 

-1 

1.168 


-2474 

10,071 

5 

7 

-2 

1.280 


-4948 

11,339 

5 

6 

-4 

1.441 


-12369 

13,809 

6 

5 

-10 

1.755 

2 

7042 

0 

1 

4 

1 

0 


' 0 

15,157 

5 

5 

0 

1 


-1761 

15,665 

6 

5 

-.25 

1.034 


-3521 

16,096 

6 

5 

-.50 

1.062 


-7042 

16,958 

6 

5 

-1 

1.119 


-14085 

18,245 

7 

4 

-2 

1.204 


-28170 

19,866 

7 

4 

-4 

1.311 


-70424 

22,989 

9 

3 

-10 

1.517 

• 

3 

247 

0 

1 

8 

1 

0 


0 

4089 

1 

6 

0 

1 


-62 

4569 

2 

8 

-.25 

1.118 


-124 

4917 

2 

8 

-.50 

1.203 


-247 

5427 

3 

8 

-1 

1.327 


-494 

6046 

3 

8 

-2 

1.479 


-988 

7033 

4 

8 

-4 

1.720 


-2470 

8868 

4 

7 

-10 

2.196 
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TABLE II 
(Continued) 


Calculated Data for Interaction 
Example Configurations - Axial 
Compression and Internal Radial Pressure 


CONFIGU- 

RATION 

N 

y 

N 

X 

m 

n 

N 

I 

N 

y o 

N 

X 

N 

X 

o 

4 

1065 

0 

1 

6 

1 

0 


0 

8347 

3 

7 

0 

1 


- 266 

8661 

4 


-.25 

1.038 


- 533 

8948 

4 


-.50 

1.072 


-1065 

9523 

4 


-1 

1.141 


-2131 

10,265 

5 

1 

-2 

1.230 


-4262 

11,303 

6 

6 

-4 

1.354 


-10654 

13,117 

7 

5 

-10 

1.571 

5 

250 

0 

1 

8 

1 

0 


0 

4301 

1 

6 

0 

1 


- '62 

4938 

2 

8 

-.25 

1.148 


- 125 

5290 

2 

8 

-.50 

1.230 


- 250 

5848 

3 

8 

-1 

1.360 


- 500 

6474 

3 

8 

-2 

1.505 


- 1000 

7360 

4 

8 

-4 

1.711 


- 2499 

8947 

5 

7 

-10 

2.080 
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Figure 3 - Stability Interaction Results for 
Combined Axial Compression and 
Internal Radial Pressure 
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SECTION 4 
DESIGN CURVES 

4.1 FIRST APPROXIMATIONS 

From the contents of the foregoing sections, it can be readily ap- 
preciated that rather complex computations are required before one can 
plot accurate interaction curves for eccentrically stiffened circular 
cylinders. Therefore, in the early stages of design, it might often be 
helpful to make crude estimates of the true interaction behavior. For 
cases of 

(a) Axial compression plus radial crushing pressure 

(b) Axial compression plus overall bending 

(c) Axial compression plus in-surface shear (N ) 

_ xy 

(d) Overall bending plus in-surface shear (N^^) 

the simple straight-line plot of Figure 4 could be used for such purposes. 

The parametric studies of SECTION 3 and reference 8 tend to indicate that 
this crude approximation would give conservative estimates for combinations 
(a) and (b), respectively. In addition, the limited test data of reference 
13 indicate similar conservatism for combination (d). 

4.2 IMPROVED APPROXIMATIONS 

In order to obtain improved accuracy over that afforded by Figure 4, 
one could make use of standard dimensional analysis concepts to arrive at 
approximate interaction curves for eccentrically stiffened cylinders. To 
understand the approach which might be taken in this connection, a discussion 
is given here for the particular case of axial compression acting in com- 
bination with radial crushing pressure. For this case, one can apply the 
theoretical solution of Block et al. [l]. The major difficulty encountered 
in attempting to develop related interaction curves arises out of the large 
number of independent variables required for the description of particular 
configurations. From equations (2-10) and (2-ll), it can be seen that, 
for cylinders made of a single isotropic material, the interaction behavior 
is a function of the following: 
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Figure 4 - First-Approximation Interaction Design Curve 
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N 

x 

N 

y 

E 

G 

V 


(A/d) 

(A r /*> 

M 

M 

(J /d) 
s 

(J /X) 
r 


R 

a 

t 

z 

s 

z 

r 


As a first step toward practical simplification of the inherent difficulties, 
it is reasonable to take p = .30 for all cases. One can then employ the 
relationship 


G 


E 

2 ( 1 +p ) 


(4-1) 


to express G in terms of E. Furthermore, for many structures, the 
torsional constants J g and will not play very significant roles. There- 
fore, in the interest of simplicity, one might choose to set both of these 
values equal to zero. In fact, this is probably a prudent practice even 
where simplification is not an objective. 

As a result of the several foregoing considerations, the array of 
variables involved in the interaction analysis reduces to the following: 


N 

x 

N 

y 

E 

R 



a 


t 

z 


z 


r 


For added convenience, 
rewritten as 


note that 


the ratios 
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where , 

I = Centroidal moment of inertia of single stringer (no 

S 

basic cylindrical skin included). 

I r = Centroidal moment of inertia of single ring (no basic 
cylindrical skin included). 

Because of these relationships, the above list of relevant variables can 
be revised to the following: 


N 

x 

N 

y 

E 

R 


(A /d) 
s 

a 

(A //) 
r 

t 

(I /d) 

s 

z 

(I /D 
r 

2 

] 


At this point it becomes helpful to apply the Buckingham Pi Theorem [l4] 
which constitutes the primary dimensional analysis concept of interest to 
the present discussion. For this purpose it is first noted that the above 
reduced listing includes 12 variables which only involve the two basic 
dimensions of force and length (pounds and inches, for example). From the 
Buckingham Pi Theorem,it therefore follows that the interaction behavior 
can be expressed in terms of 10 ( = 12-2) dimensionless ratios which may be 
chosen as 
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Thus the critical value 


for (N /ER) 

y 


can he 


expressed as follows: 



A A I 
a _s _r s 

R * dR * iR *' n 3 ’ 

dK 





(4-3) 


One 


can then proceed by introducing the quantities 



and N 


y 


o 


where 


= Critical value of axial compressive loading when 
° acting alone. 

= Critical value of circumferential compressive loading 
° when acting alone. 

Simple algebraic operations then lead to the following result: 



(4-4) 
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Note that the number of relevant variables is still unwieldy for the 
plotting of design curves. However, for selected values of the ratios 



it should be possible to select reasonable, practical magnitudes for the 
remaining ratios. Approximate interaction curves could then be plotted 
which would provide more realistic preliminary estimates than can be 
obtained from Figure 4. Such procedures were followed in the development 
of the curves given in Figure 5. These plots are furnished here only as 
examples of what might be accomplished along these lines. Within the scope 
of the investigation covered by this volume, it did not prove possible to 
refine this approach to a sufficient degree to obtain practical curves for 
direct application to actual structures. Some of the values selected in 
the generation of these curves resulted in rather unrealistic situations. 
Nevertheless, the plots do demonstrate an approach which might be further 
developed in the future. Such development should include further study 
to arrive at reasonable practical ratio values. In addition, additional 
consideration might be given to alternative formats for the data presentation 
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R/t = 500 



1 * 
-±r = lxlO* 6 
dR 3 


Stringers Outside, Frames Inside 



O .2 .4 .6 .8 1*0 


N 



Figure 5(a) - Sample Improved-Approximation 
Interaction Design Curves for 
Combined Axial Compression and 
External Radial Pressure 
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Figure 5(b) - Sample Improved - Approximation 
Interaction Design Curves for 
Combined Axial Compression and 
External Radial Pressure 


4-8 

GENERAL DYNAMICS CONVAIR DIVISION 


GDC-DDG-6 7-006 


R/t = 1000 



i = ,4X10" 6 
dR 3 


Stringers Outside, Frames inside 


R 

y 




Figure 5(c) - Sample Improved- Approximation 
Interaction Design Curves for 
Combined Axial Compression and 
External Radial Pressure 
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R/t = 1000 | = 2.0 

H 


— I 2xlCf 6 
dR 3 


Stringers Outside, Frames Inside 



0 .2 .4 .6 .8 1.0 



Figure 5(d) - Sample Improved-Approximation 
Interaction Design Curves for 
Combined Axial Compression and 
External Radial Pressure 
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R/t = 5000 | = 2.0 

K 


-ir = 1.5xl0" 8 
dR 3 


Stringers Outside, Frames Inside 




Figure 5(e) - Sample Improved- Approximation 
Interaction Design Curves for 
Combined Axial Compression and 
External Radial Pressure 
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1 

R/t = 5000 f- = 2.0 — % = 7.5x10 

R dR 




Figure 5(f) - Sample Improved-Approximation 
Interaction Design Curves for 
Combined Axial Compression and 
External Radial Pressure 
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SECTION 5 

DIGITAL COMPUTER PROGRAM 

This section presents the essential features of General Dynamics 
Convair digital computer program numbered 3962. This program was de- 
veloped to implement the theoretical solution developed by Block, et al. 

Cl] for the buckling of eccentrically stiffened circular cylinders sub- 
jected to axial load and/or radial pressure. The program allows for the 

input of N or N (+ = compression, - = tension) along with the m and n 
y x 

values to be screened. 

m = Number of axial half-waves in buckle pattern, 
n = Number of circumferential full-waves in 
buckle pattern. 

The (or N ) value is computed for each ra-n combination specified by 

the user. The lowermost computed N (or N ) is always retained throughout 

x y 

the process and the final overall minimum value is printed out as the critical 
loading. Judgement must be exercised by the user to insure that the actual 
minimum- strength buckle pattern does not lie outside the range of m and n 
values screened. The input format is shown in Figure 6. Symbols are listed 
in Table III. A detailed, card-by-card description of the input follows 
below. Runs may be stacked. 


CARD TYPE 1: One card per run. 

Enter PROBLEM IDENTIFICATION (any alphanumeric 
characters) in columns 1-72. 

CARD TYPE 2: One card per run. 

Enter NCASES (number of cases) as right adjusted integer 
in columns 1-5. 

Enter PRNTOP (printout option) as right adjusted integer 

in columns 6-10. If PRNTOP = 0 (or blank), only output 

corresponding to minimum N are printed. If PRNTOP = 1 

(orqfcO), all output associated with each m and n 

combination screened are printed including the A. . values 

— J 

used in equation (2-10). 
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Figure 6 - Input Format - Program 3962 
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CARD TYPE 3: One card per case* 

Enter CASENO (case number) as right adjusted alphanumeric 
characters in columns 1-5. 

Enter LOADIN (loading condition) in columns 6-10 with 
the symbols: 

COMBY if N value input on CARD TYPE 3. 

_y 

COMBX if N x value input on CARD TYPE 3. 

Enter MINM (minimum m value for screening) in columns 
11-20 (E10.5). MINM > 0. 

Enter MAXM (maximum m value for screening) in 
columns 21-30 (£10.5). 

Enter DELTAM (Am) in columns 31-40 (E10.5). 

Enter MINN (minimum n value for screening) in 
columns 41-50 (E10.5). MINN > 0 for COMBX. 

MINN Si 0 for COMBY. 

Enter MAXN (maximum n value for screening) in 
columns 51-60 (E10.5). 

Enter DELTAN (An) in columns 61-70 (E10.5). 

Enter NY (N ) or NX (N ) in columns 71-80 (E10.5). 

y x 

When NY is input, the program computes the critical 
NX value. When NX is input, the program computes the 
critical NY value. (NX and NY are positive for com- 
pression and negative for tension) . 

CARD TYPE 4: One card per case. 

Enter RADIUS (R, cylinder radius measured to mid- 
surface of skin, in.) in columns 1-10 (E10.5). 

Enter LENGTH (a, overall length of cylinder, in) 
in columns 11-20 (E10.5). 
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Enter STRSP (d, stringer spacing, in.) in columns 
21-30 (E10.5). 

Enter FRAMSP (f, frame spacing, in.) in columns 
31-40 (E10.5). 

CARD TYPE 5: One card per case. 

Enter EX (E , lbs/in) in columns 1-10 (E10.5). 
x’ 

Enter EY (E , lbs/in) in columns 11-20 (E10.5). 

y 

Enter DX (D , lb-in) in columns 21-30 (E10.5). 
x J 

Enter DY (D , lb-in) in columns 31-40 (E10.5), 

y 

Enter GXY (G , lbs/in) in columns 41-50 (E10.5). 
xy’ 

Enter DXY (D , lb-in) in columns 51-60 (E10.5). 
xy’ 

CARD TYPE 6: One card per case. 

Enter PRIMMX (p ') in columns 1-10 (E10.5). 
r x 

Enter PRIM4Y (p •) in columns 11-20 (E10.5). 

y 

Enter MX (p ) in columns 21-30 (E10.5). 
x 

Enter MY (p ) in columns 31-40 (E10.5). 

y 

CARD TYPE 7: One card per case. 

Enter ES (E , psi) in columns 1-10 (E10.5). 

s 

Enter GS (G , psi) in columns 11-20 (E10.5). 
s 

Enter AS (A , in 2 ) in columns 21-30 (E10.5). 
s 

Enter IOS (I , in 4 ) in columns 31-40 (E10.5). 

°s 4 

Enter JS (J , in ) in columns 41-50 (E10.5). 

S 

Enter ZBARS (z , in) in columns 51-60 (E10.5); positive 
s 

for external stringers; negative for internal stringers. 
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CARD TYPE 8: One card per case. 

Enter ER (E^, psi) in columns 1-10 (E10.5)„ 

Enter GR (G , psi) in columns 11-20 (E10.5). 
o 

Enter AR (A . in ) in columns 21-30 (E10.5). 
r 

Enter IOR (T , in^) in columns 31-40 (E10.5). 
o 

r . 

Enter JR (J , in ) in columns 41-50 (E10.5). 
r 

Enter ZBARR (z^, in) in columns 51-60 (E10.5); positive 
for external rings; negative for internal rings. 

A sample input coding form is shown in Figure 7. A sample output listing 
is given in Figure 8. A basic flow diagram is presented in Figure 9 and 
a Fortran listing of the program is shown in Table IV. 
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TABLE III - Program 3962 Notation 


PROGRAM 

REPORT 


NOTATION 

NOTATION 

DESCRIPTION 

AR 

A 

Cross-sectional area of circumferential 

r 

stiffener, in . 

AS 

A 

Cross-sectional area of longitudinal 


s 

2 

stiffener, in . 

CASENO 

- 

Case number. 

COMBX 

— 

Indicates N^ value input, solve for 

N . 

y 

COMBY 


Indicates N value input, solve for 

y 

N . 

X 

D 

- 

See STRSP. 

DX 

D 

X 

Bending stiffness of skin in longitudinal 
direction, lb-in. 

DY 

D 

y 

Bending stiffness of skin in circum- 
ferential direction, lb-in. 

DXY 

D 

xy 

Twisting stiffness of skin, lb-in. 

DELTAM 

Am 

Value by which m is incremented. 

DELTAN 

An 

Value by which n is incremented. 

ER 

E 

r 

Young's modulus for circumferential 
stiffener, psi. 

ES 

E 

s 

Young's modulus for longitudinal 
stiffener, psi. 

EX 

E 

X 

Extensional stiffness of skin in 
longitudinal direction, lbs/in. 

EY 

E 

y 

Extensional stiffness of skin in 
circumferential direction, lbs/in. 
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TABLE III 

- Program 3962 Notation 

PROGRAM 

REPORT 

(Continued) 

NOTATION 

ERARL 

NOTATION 

EA/I 
r r 

DESCRIPTION 

FRAMSP 

FMPIA 

FMXMY 

GMXMY 

a 

nm/a 
1-p ’p ' 

x 

1-p V 
x y 

Frame spacing, in. 

GR 

G 

r 

Shear modulus for circumferential 
stiffener, psi. 

GS 

G 

s 

Shear modulus for longitudinal 
stiffener, psi. 

GXY 

G 

xy 

In-plane shear stiffness of 
skin panel, lbs/in. 

IOR 

I 

o 

r 

Moment of inertia of circumferential 
stiffener cross-section about middle 
surface of skin, in 4 . 

IOS 

I 

0 

s 

Moment of inertia of longitudinal 

stiffener cross-section about middle 

4 

surface of skin, in . 

JR 

J 

r 

Torsional constant for circumferential 
stiffener, in 4 . 

JS 

J 

s 

Torsional constant for longitudinal 
4 

stiffener, in , 

LENGTH 

a 

Length of stiffened cylinder, in. 

LOADIN 

- 

Option describing input loading: 
COMBY or COMBX. 

MAXM 


Maximum value of m used. 
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TABLE III - Program 3962 Notation 
(Continued) 


PROGRAM 

REPORT 


NOTATION 

NOTATION 

DESCRIPTION 

MINM 

- 

Minimum value of m used. 

MAXN 

- 

Maximum value of n used. 

MINN 

- 

Minimum value of n used. 

MX 

V 

X 

Poisson's ratio for bending of skin 
in longitudinal direction. 

MY 

V 

Poisson's ratio for bending of skin 
in circumferential direction. 

NBAR 

N 

N or N (+ in compression) Ibs/in. 
x y 

NCASES 

- 

Number of cases. 

NOP 

- 

NOP = 1 used for COMBY 
NOP = 2 used for COMBX 

NXORNY 

N or N 
x y 

N or N input, Ibs/in. 
x y 

PRIMMX 

*x 

Poisson's ratio for extension of skin 
in longitudinal direction. 

PRIMMY 

* 7 ' 

Poisson’s ratio for extension of skin 
in circumferential direction. 

PRNTOP 


Printout option. If PRINTOP = 0 (or 
blank) , prints only output associated 
with minimum value of N calculated. If 
PRNTOP =£0, prints output for all com- 
binations of m and n values screened 

PROBID 

- 

Problem identification. 

RADIUS 

R 

Radius of cylinder, measured to mid- 
surface of skin, in. 
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TABLE III - Program 3962 Notation 
(Continued) 


PROGRAM 

REPORT 


NOTATION 

NOTATION 

DESCRIPTION 

STOREM 

- 

m for minimum N . 

STOREN 

- 

n for minimum N . 

STRSP 

d 

Stringer spacing* in. 

TEMP 

- 

Relative minimum N . 

ZBARR 

z 

r 

Distance from middle surface of skin to 
centroid of circumferential stiffener, 
in., positive if stiffener is outside. 

ZBARS 

z 

s 

Distance from middle surface of skin to 
centroid of longitudinal stiffener in., 
positive if stiffener is outside. 
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FORTRAN CODING AND DATA FORM 



1 -MMJEH 


m 

Si 


v 

i**? L 

EWO 



c< 

BE* 

i-.jS 


’X 

'Z 





■ 

bbhihhhhhh 

■ 

OAT 


1m$F .j 

* 

W 

/ 



1 

2 


r < 

i* 

Q 

Q 

0 

D 

B 

□ 

C 

Q 

B 

ft 

D 

i 

I 

B 

c 

B 

B 

B 

B 

B 

B 

B 

B 


C 

B 

E 


E 

B 

E 

B 

3 

3 

3 

i 

3 

3 

3 

3 

S 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

31 

B 

1 

31 

B 

B! 


B! 

3 

B! 

B 

3 


i 



1 

1 

1 

1 


1 

1 

5 

D 

ft 

5 

ft 

ft 

1 



ft 

9 


5 

ft 




1 







1 


1 

1 



1 


1 

1 




1 

1 

1 




1 



1 








1 

1 

1 








1 

1 

1 





l 

2 



1 

1 

ft 

5 


7 

\ 

ft 



1 

E 


E 


1 

1 

1 



& 

5 


E 

1 







B 


E 

1 



1 

i 


1 



1 

1 


1 



1 


5 


E 








n 


E 


1 



1 



Si 




i 



2 

3 



l 

\ 

8 

l 


1 

1 

1 



ft 

ft 


ft 


1 

1 

1 




5 


r, 

E 







ft 


S 

1 






1 



1 

1 

1 





1 



1 








a 

1 






1 










3 

4 



1 

E 

8 

5 


r 

tj 

E 

ft 



1 

ft 


ft 


r 

rtf 

ft 

ft 




ft 


E 

1 







ft 


ft 

1 





1 




E 

1 

5 

E 



E 

ft 

5 


ft 








i 

1 






1 










4 

5 



1 

1 

8 



a 

1 

1 



1 

1 


ft 


a 

1 

1 




1 


E 

5 









E3 

5 





1 




1 

1 

1 




1 


1 


1 















1 










5 

6 



Z 

ft 

i 

ft 


s 

ft 

ft 



E 

ft 


ft 


s 

ft 







ft 

ft 









ft 

ft 





1 




ft 

ft 

ft 




8 




1 








a 

1 
















6 

7 



r 


8 

ft 


P! 

rtf 

g 

ft 



E 

ft 


ft 


E 

ft 

ft 




1 


ft 

ft 









ft 






1 




ft 

ft 

ft 




8 


ft 


1 








a 
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7 

8 
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1 

1 


1 

1 
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1 





1 

1 

1 




1 



1 















1 




1 

1 

1 




1 


1 


1 








a 

1 






1 










8 

9 



i 


1 

1 


1 


1 



1 






1 

1 




1 


















1 




1 

1 

1 




1 


1 


1 








a 

1 






1 










9 

10 



I 


1 

1 


1 


1 



1 






1 

1 






















1 




1 

1 

1 




1 


1 


1 








a 

1 






1 










10 

n 



i 


1 

1 


1 


1 



1 






1 
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1 

1 
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1 
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CASE NO. 

MIN M 

MAX M 

DELTA M 

MIN N 

MAX N 

DELTA N 

1 

1.00GC1E m 

2.5GOCE 01 

l.OOOOE 00 

-n. 

2.5GPDE 01 

l.OOOCE DO 


CYLINDER 

RADIUS 

OVERALL 

LENGTH 

STRINGER 

SPACING 

FRAME 

SPACING 


N BAR Y 


3.8600E 01 

7.2000E 01 

2.480OE 00 

6.0C00E CO 


2.0000E 00 


E SUB X 

E SUB Y 

D SUB X 

D SUB Y 

G SUB XY 

D SUB XY 


1.5000E 06 

2.0000E 06 

2.5000E 02 

3. 0000 E 02 

2.0000 E 05 

2.0CC0E 02 


MU PRIME 
SUB X 

MU PRIME 
SUB Y 

MU SUB X 

MU SUB Y 




2.5000 £-01 

3.5000E-01 

2. 50Q0E-C1 

3. 5000E-01 


Z BAR 


E SUB S 

G SUB S 

A SUB S 

I SUB OS 

J SUB S 

SUB S 

01 

1 

M 

3. 0000 E 07 

1.2000E 07 

2.000OE-G2 

5.noC0E-C3 

4.C000E-03 

5.0000E-C2 

M 







- 

E SUB R 

G SUB R 

A SUB R 

I SUB OR 

J SUB R 

Z BAR 
SUB R 


2. 5000 E 07 

l.OOOCE 07 

4.CQOOE-OZ 

1.0000E-02 

6.0000E-03 

-1.0OCQE-O1 



CRITICAL 

COMBINED LOAD VALUES 


NUMBER OF 

NUMBER OF 

AXIAL LOAD 

HOOP LOAD 

LONGITUDINAL 

CIRCUMFERENTIAL 

BAR SUB X 

BAR SUB Y 

HALF WAVES 

HALF WAVES 

LBS PER IN 

LBS PER IN 

4.0000E 00 

1.4000E 01 

7.8588E 03 

2.QOO0E 00 

3 .0000 E 00 

1.40COE 01 

8.0546E 03 

2.P00CE 00 

5.00C0E 00 

1 .4000 E 01 

8.3383E 03 

2.PU00E 00 

4.iH)00E 00 

1.2COOE 01 

8.1912E 03 

2 • 00(50 E PC 

4.0000E 00 

1.6C0OE Cl 

8.3470E 03 

2.0000E 00 


Figure 8 - Sample Output Listing - Program 3962 
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Figure 9 - Flow Diagram - Program 3962 
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TABLE IV - Fortran Listing - Program 3962 


SIbFTC MAIN 
COMMON 

1 

2 

3 

DIMENSION 


LIST 

/INPT/ 


CaSEnO»LOADIN»mINM' ma XM»DELTAM»MiNN»MaxN»OELTaN, 

RAOIuS»LENGTH»STRsP»FRAMSPf EX , EY , DX * QY » GXY » DXY , 
PR I MMX . PR IMM Y , MX » M Y » ES . GS * AS , IOS , JS . ZBaRS t ER . GR » 
AR » I OR » JR » ZBARR , NoP > P I * PRNTOP , NXORMY 
PR0BID<12> 


REAL MlNMfMAXMfMiNN»MAXN,LENGTH»Mx»MY , IOS. JS. I OR , JR * NRaR » NXORNY 

integer combx, comby 

DATA COMBX . COMBY /5HC0MBX » 5HC0MB Y/ 

DATA Pi /3.14IB926/ 

100 READ (5.101) PROBID 

101 FORMAT (12A6> 

WHITE 16.201) PROBID 

201 FORMAT (lHl»20(/> .58X. 16HLANGLEY SOLUTION 1 0 ( / ) . 3 OX » 12a& ) 

READ ( b» 301 ) NCASES. PRNTOP 
301 FORMAT (IB.F5.0) 

DO 1000 N=1.NCASLS 

call input 

IF (PRNTOP. EQ. 0. ) GO TO 310 
WRITE (6.305) CASENO 

305 FORMAT ( 1H1 » 50X * 2sH INTERMEDIATE VALUES FOR CASE.aS) 

310 NOP=0 

IF (LOADIN.EQ. COMBY) NOP=l 
IF (LOADIN.EQ. COMBX) N0P=2 
IF (NOP.NE.O) GO TO 400 
Vv R 1 TE (6.350) CASENO 

350 FORMAT (/// 87X » 3lHOPT I ON GIVEN INCORRECTLY. CASE A3. gH DELETED) 
CO TO 1000 
400 CALL EQUATN 

500 call final 

t>00 CALL OUTPUT 
1000 CONTINUE 
GO TO 100 
END 

SIbFTC NDAK LIST 
REAL FUNCTION 

1 


N8AR(EX.FMxMY.ES. AS.D.FMPIA.GXY.FN r.PRImmY.R.ZBARS. 
EY.ERaRL.ZBARR.DX.I0S.DXY.6S. JS.GR. jR.L.DY.ER. 
IOR.N.NOP.M.PRnTOP.'NXORNY.mY.gMXMY) 

Real IOS. JS# JR »L*iOR»N»M. NXORNY 
real my 

ah = ( Ea / FMXMY + ES*AS/D) * FMpIA**2 + GXY*FNR**2 
CPR1MMY*EX/F ; MXMY + GXY) * FmPIA * FNR 

(l./R)*(PRiM:MY*EX/FMXMY)*FMplA + ES*AS/D*ZbARS*FMPI A**3 


A12 = 
A13 = 
A22 - 
A23 = 
A33 = 


GXY *FMP I A**2 + (EY/FMXMY + £ RARE ) * FNR**2 
(l./R) * (LY/FMXMY + ERARL) * FNR + ERARL*ZBARR*FNR**3 
(DX/GMXMY + ES*I0S/D)*FMPIA**4 + ( 2 • *MY*DX/GMXMY 

1 + 2 • *DX Y + Gb* JS/D + GR*JR/L } * FMPIA**2 * FNR**2 

2 + (DY/GMXMY + £R*I0R/L> *FNK**4 f (1»/R**2)* (EY/FMXMY 

3 + ERARL) + (2.*ERARL*ZBARR*(N**2/R**3) ) 

GO TO (1000.2000), NOP 
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TABLE IV - Fortran Listing - Program 3962 
(Continued) 


C N0P=1 

1000 DIV = FMPIA**2 

SUB = FNR**2 * NXORNY 
GO TO 4000 
C 

C NOP=2 

2000 DIV = FNR**2 

SUB = FMPIA**2 * NXORNY 
C 

4000 NBAR = (A33 + ( ( ( A12*A23) - ( A13*A22 > ) / ( ( A11+A22 ) -A12**2 ) ) *A13 

1 + .( ( ( (A12*Al3)-(All*A23J )/( (A11*A22)-A12**2) )*A23 ) - SUB) 

2 / DIV 

IF (PRNTOP.EQ.O.) RETURN 

WRITE (6 * 6000 ) M»N » NBAR » All » Al2 . Aj.3» A22 ♦ A23. A33 
6000 FORMAT (/// 11X»5 hM = » lPEU.3.2aX»5HN = » lPEll.3*28X»6HN BAR= . 

1 E16.9 // HXfbHAll =,1PE17.9.22X»5HA12 =.lPEl7.9f 22X. 

2 5HA13 =»1PE17.9» // 11X.5HA22 =. 1PE17.9.22X»5HA23 =.1PE17.9, 

3 22X . 5HA33 =»1PE17.9) 

RETURN 

END 

5IBFTC INPUT LIST 

SUBROUTINE INPUT 

COMMON /INPT/ CaSEnO* LOADING MINm»MAXM,DELTAM»MiNN»MaxN»DELTAN» 

1 radius » length . strsp * fr amsp > ex . ey , dx . dy . gx y . dx y . 

2 PR I MMX > PR I MM Y » MX ' MY » ES . GS . AS . I OS , JS . ZB ARS . ER . GR » 

3 AR# I0R» JR*ZBARR»NoP»PI >PRNTOP, NXORNY 

REAL M INM » MAXM # MlNN»MAXN» LENGTH rMx» MY » IOS> JS» IOR» JR NXORNY 

C 

C SUBROUTINE TO READ CARD INPUT FOR EACH CASE. 

READ (5.101) CA5EN0»L0ADIN,MlNM>M A XMfDELTAM»MINN,MAXN»DELTAN 
* » NXORNY 

ioi format (2A5*7eio.5) 

read (5.201) RADluSfLENGTH»STRSP»FPAMSP 
201 FORMAT (4E10.5) 

READ (5.301) EX»£y»DX.DY.GXY*DXY 
301 FORMAT (bElO.S) 

READ (5.201) PRIMMX.pRIMMY.MX.MY 
READ ( 5» 301 ) ES»GS»AS» IOS* JSrZBARs 
READ (5.301) ER.GR.AR.IOR. JR.ZBARR 
RETURN 
END 

■fclBFTC EGUATN LIST 

SUBROUTINE EQUATN 

common /inpt/ caseno*loadin»miimm»maxm.deltam»minn.maxn»deltan» 

1 Radius » length. strsp* framsp. ex. ey.dx.dy.gxy.dxy. 

2 PrIMMX>PRIMMY»MX'MY»ES»GS»AS» I.OS, JS.ZRARS.ER.GR. 

3 AR» I0R» JR»ZBARR»NoP'PI .prntop. nxorny 
COMMON /INTER/ FMxMY . ERARL . TEMP . S jOREM » STOREN . GMxMY 
COMMON /RENAME/ A.D.L.R 

REAL M I NM . MAXM . Ml NN . MAXN . LENGTH . Mx . MY , IOS . JS » IOR , JR . NBAR 
REAL L*M.N.NBARX.nBARY»NXORnY 
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TABLE IV - Fortran Listing - Program 3962 
(Continued) 


R=RAUIUS 

A=LEN6TH 

D=STRSP 

L=FKAMSP 

FMXMY = i.-PRlMMX*PRIMMY 
GMXMY = l.-MX+MY 
ERARL = LR+AR/L 
1CNT=0 

NsMINN-DELTAN 
ioo n=n+uEltan 

IF IN.GT.MAXN) GO TO 6000 
F NR = N/R 
M=MiNM-DELTAM 
200 M-M+DELT AM 
1CNT=1CNT+1 

IF tM.GT.MAXM) GO TO 5000 

fmpia = m*pi/a 

60 TO ( 1000 » 2000 ) » NOP 

c 

c NOP=i SOLVE FOR N BAR X 

1U00 IMBARX - NBARIEX, FMXMY, ES, AS, D»FmPiA»GXY,FNR»PRIMmY»R»ZBARS,EY, * 

1 ERARl , ZtJ ARR.DX. IOS,DxY»GS» JS»6R» JR.L'DY.Fr. I OR, N, NOP, 

2 MtPKNTOPfNXORNY»MYfGMXMY) 

IF (ICNT.NE.l) GO TO 1100 

tlmpsnbarx 

ST OKEM=M 
STOREN=N 
GO TO 4000 

1100 IF ( NBARX • GE • TEMP ) GO TO 4000 
1 £MP— NBARX 
STOREM=M 
ST OKEN = N 
60 10 4000 

NOP=2 SOLVE FOR N BAR Y 
2000 NBARY = NBAR(EX.FXXMY.ESf AS. D. FMPiA. 6 xY,FNR.PRIMmY»R»ZbARS»EY. 

1 ERARl.2BARR.OX.IOS.OxY.GS.JS.GR. JR, l.DY.ER. ior.n. nop , 

2 M.PRnTOP.NXORNY. MY, GmXMY) 

IF (ICNT.NE.l) GO TO 2100 
TEMP=NBAKY 

STORtM=M 
STORLN=N 
60 10 4000 

2100 IF ( NBARY. GE. TEMP) GO TO 4000 
T EMP— NBARY 
STOREM=M 
ST0REN=N 

4000 60 TO 200 
5000 60 TO 100 
6000 RETURN 
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TABLE IV - Fortran Listing - Program 3962 
(Continued) 


END 

SIBFTC FINAL LIST 

SUBROUTINE FINAL 
COMMON /INPT/ 


CASENO»LOADlNfMINM»MAXM,DELTAM»MiNN»MAxNrDELTAN# 

1 RADIUSfLENGTHfSTRSP»FRAMSP»EX»EY,DXrDY»GXY»DXY» 

2 PRlMMXrPRlMMYrMX'MY'ESrGSr ASrIOSr JSrZBARSrERrGRr 

3 Ar# ior» jr»zbarr»nop»pi r prntop, nxorny 

COMMON /OUT/ X { 5) > Y ( 5) r VALUE ( 5 ) 

COMMON /INTER/ FMxMY r ERARLr TEMP r S t0REM» STOREN* GMxMY 
COMMON /RENAME/ A,DrLrR 
DIMENSION FN(S) 

REAL MlNMrMAXM»MlNNrMAXNr LENGTH »M x^MY,IOS^JS»IOR,OR»NBaR»M»N 
1 r NXORNY 

valueu)=temp 

X(1)=ST0REM 

Y(1)=2.0*STOREN 

FN(1)=ST0REN 

X(2)=X(1)-1.0 

X(3)=XU)+1. 

X(4)=X(1) 

X(5)=X(1) 

Y(2)=Y(1) 

FN(2)=FN(1) 

Y(3)=Yll) 

FN(3)=FN(1) 

Y(4)=Y(l)-2.0 
FN(4)=FN(1)-1.0 
Y(5)=Y(l)+2.0 
FN(5)=FN(l)+I.O 
DO bOO I=2»b 
FMP1A= X(I)*PI/A 
FNR = FNII)/R 
N = FN(I) 

M-X ( I ) 

VALUE ( I ) 


1 

2 

bOO CONTINUE 
RETURN 

end 

SIBFTC OUTPUT 


NBAR(Ex*FMXMY,ES*ASrD*FMPlArGXY,FNR»PRlMMY,R»ZBARS*EY 
ErARL r ZBARR r DX r I0f> r DXY r GS » JS r GR » JR r L r D Y r ER r I OR r N 
NQP r M » PRNTOP r NXORnY * MY r GMXMY ) 


LIST 


SUBROUTINE OUTPUT 

C OMMON /INPT/ C ASENO t LOAD I N r M I N m » M AXM r DELT AM r M I NN r MAxN r DELT AN r 

1 RADIUS »LENGTH»STRSP'FRAMSP» EX *EY,DXrDY»GXY,DXY» 

2 PRlMMX»PRlMMY>MX»MY»ES»GS»AS»lOSf JS»ZBARSfER*GR» 

3 AR r I OR » JR r ZBARR r NqP » P I » PRNTOP f NXORNY 

COMMON /OUT/ X ( b) r Y ( 5) r VALUE ( 5 ) 

REAL MlNMrMAXM»MlNN»MAXN»LENGTH»Mx»MY» IOSr JSr IOR r JR » NBaR r NXORNY 
WRITE <6rlOO> CASENO»MlNM»MAXMfDELTAM»MINNrMAXNrDELTAN 
100 FORMAT (1H1*1X»8HCASE N0.»9X»5HMlN M»16X»5HMAX M r 14X r 7HDELTA Mr 
1 14X r 5HMIN N» 15X * 5HMAX Nr 14x » 7HDELTA N // IX r A6 r 1P6E20 .4 ) 
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TABLE IV - Fortran Listing - Program 3962 
(Continued) 


GO TO (200,250). NOP 

200 w R i T L (6.201) RADIUS. LENGTH . STRSP , FRAMSP » NXORnY 

201 FORMAT (/// l8X#6HCYLlNDER#l2X»7HoV£RALL»13Xr8HSTRlNGER»13X» 

1 5HFRAME / l9X .6HRADIUS. 14X , 6HLENGTH. 14X» 7 hSPACInG» 12X . 

2 7HSPACING.34X.7HN BAR Y // 7X . 1P4E20 .4 , 1PE40 .4 ) 

00 TO 299 

2b0 WRITE (6.251) RADIUS. LENGTH. STRSP, FRAMSP. NXORnY 

251 FORMAT (/// 18X . ohCYLINDER . 12X . 7HqVEKALL . 13X . 8HSTRINGER » 13X . 

1 5HFRAME / 19X . 6HRADIUS . 14X , 6HLENGTH. 14X » 7hSPACInG » 12X . 

2 7HSPACING.34X.7HN BAR X // 7X . 1P4E20 . 4 . 1PE40 . 4 ) 

299 WRITE (6.300) EX * E Y . OX . DY . GXY , DX Y 

300 format (/// ibx.'/hE sub x»i3X.7he sub y,i3x.7hd sub x,i3X. 

1 7hD SUB Y.13X.8HG SUB XY.I2X.8HD SUB XY // 7X»lp6E20.4 ) 

WRUE (6.400) PR1MMX.PRIMMY.MX.MY 
400 FORMAT (/// I8X.8 HMU PRIME. 12X»8HmU PRIME / 19X»5HSUB X.15X. 

1 5HSUB Y.14X.8HMU SUB X.12X,8HMU SUB Y // 7X.1P6E20.4 ) 

WRITE (6.500) ES» OS. AS» IOS. JS.2BARS 
500 FORMAT (/// 1 19X » 5HZ BAR / 16X.7Hf; SUB S.13X.7HG SUB S.13X. 

1 7HA SUB S.13X.8HI SUB 0S.l2X.7*HJ SUB S.l4x»5HSUR S // 

2 7X . 1P6E20 • 4 ) 

WRITE (6.600) ER»GR» AR» IOR. JR.2BARR 
oOO FORMAT (/// U9X.&HZ BAR / i8X.7Hg SUB R.13X.7HG SUB R.13X. 

1 7HA SUB RM3X.8HI SUB OR»i2*' 7 HJ SUB R.l4x*5HSUB H // 

2 7X . 1P6E20 • 4 ) 

WRITE (6.700) 

700 FORMAT (b(/) . SIX . 29HCR IT IC AL COMBINED LOAD VALUES /// 

* 12X. 9HNUivBER OF, 24X, 

1 9HNUMBER OF , 24X , 10HAX I AL L 0AU » 23X , 9HH00P LOAD / I1X, 

2 12HL0NGI JuDInAL»19X»15HCIrCUMFERENTIAL*2iX.9HbAR SUB X. 

3 24X.9HBAK SUB Y / 12X,10HhALF WAVES . 23X » 10HHALF WAVES. 

4 23X.10HLBS PER IN . 23X . 1 0H L BS PER IN // ) 

GO TO (1000.2000) » NOP 

10QO WRITE (6.1001) ( X ( 1 ) » Y ( 1 ) . VALUE ( I ) » NXORNY » 1 = 1,5) 

1001 FORMAT (1X,1PE21.4.1P3E33.4//(1X» 1PE21.4, 1P3E33.4) ) 

RET URN 

2000 WRITE (6,1001) (X( I) , Y ( I ) .NXORNY. VALUE( I ) » 1=1,5) 

3000 RETURN 
End 
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